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ABSTRACT 

Context. Several mechanisms for forming brown dwarfs have been proposed, which are today believed not to be mutually exclusive. 
Among the fundamental characteristics of brown dwarfs which are intrinsically tied to their origin, multiplicity is of particular rele- 
vance. Any successful determination of the actual mass for such objects in systems is thus worthwhile as it allows one to improve on 
the characterization of the multiplicity properties (e.g, frequency, separation, mass-ratio distribution) of sub-stellar companions. 
Aims. We attempt to improve on the characterization of the properties (orbital elements, masses) of two Doppler-detected sub-stellar 
companions to the nearby G dwarfs HD 131664 and HD 43848. 

Methods. We carry out orbital fits to the Hipparcos Intermediate Astrometric Data (IAD) for the two stars, taking advantage of the 
knowledge of the spectroscopic orbits, and solving for the two orbital elements that can be determined in principle solely by astrom- 
etry, the inclination angle i and the longitude of the ascending node Q.. A number of checks are carried out in order to assess the 
reliability of the orbital solutions thus obtained. 

Results. The best-fit solution for HD 131664 yields i = 55 ± 33 deg and CI = 22 ± 28 deg. The resulting inferred true companion 
mass is then M c = 23+ 26 Mj. For HD 43848, we find i = 12 ± 7 deg and Q. = 288 ± 22 deg, and a corresponding M c = \20 + \f Mj. 
Based on the statistical evidence from an F-test, the study of the joint confidence intervals of variation of i and fl, and the comparison 
of the derived orbital semi-major axes with a distribution of false astrometric orbits obtained for single stars observed by Hipparcos, 
the astrometric signal of the two companions to HD 131664 and HD 43848 is then considered detected in the Hipparcos IAD, with a 
level of statistical confidence not exceeding 95%. 

Conclusions. We constrain the true mass of HD 131 664b to that of a brown dwarf to within a somewhat statistically significant degree 
of confidence (~ 2— cr). For HD 43848b, a true mass in the brown dwarf regime is ruled out at the 1 - cr confidence level. The results 
are discussed in the context of the properties of the (few) close sub-stellar and massive planetary companions to nearby solar-type 
stars and their implications for proposed models of formation and structure of massive planets and brown dwarfs. 

Key words, stars : individual : HD 131 664, HD 43 848 - planetary systems - astrometry - methods : data analysis - methods : numerical 
- methods: statistical - stars: low mass, brown dwarfs - stars: statistics 



1. Introduction 

Similarly to the detection of the fi rst Jupiter-mass (M.r) p lanet 
orbiting a star other than the Sun dMavor & Quelozlll995l) . the 
first unambigu ous discovery of a brown dwarf] dates b ack ~ 14 
years (e.g.. (Nakajima et al][l99H iRebolo et al.lll995l) . Several 
mechanisms for forming brown dwarfs have been proposed: 
Turbulent fragmentation of molecular clouds, fragmentation of 
massive prestellar cores, protoplanetary disk fragmentation, dy- 
namical ejection of protostellar e mbryos, and photoeros ion of 
protostellar cores (for a review, see lWhitworth et al. 2007). Such 
mechanisms are today believed not to be mutually exclusive, and 
they all likely operate in nature. To determine their relative con- 
tribution to the overall brown dwarf population is a major the- 
oretical challenge. To date, agreement is still to be reached on 
fundamental issues such as the minimum mass for brown dwarfs, 
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1 We operationally adopt the commonly used definition of brown 
dwarf as a deuterium-burning sub-stellar object with mass in the ap- 
proximate range 15 - 80 Mj 



and how brown dwarfs might be distinguished from planets. For 
example, if brown dwarfs are identified as objects that form 
like stars do, on dynamical timescales by gravitational insta- 
bility, regardless of the formation locus (a molecular cloud or 
a marginally unstable protoplanetary di sk) and if the minimum 
mass for core collapse is a few My (see lWhitworth et ail l2007). 
then significant overlap between the mass range of brown dwarfs 
and planets would occur, and a simple mass cutoff (such as the 
one adopted here) may not apply. Given the significant num- 
ber of open questions in this research field, it is thus essential 
for theory to accurately reproduce, and for observations to care- 
fully determine, the ensemble properties of brown dwarfs, in- 
cluding the brown dwarf initial mass function, the young brown 
dwarfs kinematics, distribution, and disk frequency, and the bi- 
nary statistics of brown dwarfs across a wide ra nge of primary 
masse s a nd orbital separations (for a review see iLuhman et al.l 
120071 and lBurgasser et alj|2007l) . 

Among the fundamental characteristics of brown dwarfs 
which are intrinsically tied to their origin, multiplicity proper- 
ties (e.g, frequency, separation, mass-ratio distribution) are of 
particular relevance, as together with the possibility of determin- 
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ing the actual internal composition and atmospheric features of 
individual objects, they constitute one of the few ways to obser- 
vationally distinguish between planets and brown dwarfs in the 
possible overlap region in mass. For example, the frequency of 
close (a < 5 AU) stellar (M c > 0.08 M ) companions to nearby 
(d < 50 pc) solar-type stars is 13 + 3% dDuquennov & Mayor! 
1199 lb . while in the same range of separations the frequency 
of gi ant planets (M c < 15 My) is known today to be about 
7% dMarcv et al.l2008l) . Brown dwarfs, on the other hand, appear 
conspicuously absent in the datasets collected by decade-long, 
high-precision radial-velocity surveys of thousands of bright 
normal stars, despite the fact that their large RV signals would 
have been easily spotted. Initial claims that the frequency of 
close brown dwarf companions seemed to b e in fair agreemen t 
with a constant distribution of mass ratios (M ayor et al.l fl992) 
where later dismissed by studies that showed, based on a combi- 
nation of radial velocity measurements and Hipparcos astromet- 
ric observations, how m ost of these putativ e sub-stellar compan- 
ions where in fact stars (Hal bwachs et a l. 2000). Indeed, among 
close companions, brown dwarfs appear outnumbered by stars 
and planets by factors of ~ 100 and ~ 50, respectively, wit h typi- 
cal frequency estimates of * 0.1% (Marcy & Butlerl2( )00). Only 
about half a dozen close companions with minimum masses in 
the brown dwarf regime are known today around bright, nearby 
solar-type stars It is conceivable that observational biases might 
somewhat contribute to a reduction of the discovery rate of 
brown dwarfs with respect to planets (massive companions be- 
ing assigned typically lower priorities than lower-mass planets 
in Doppler surveys hard-pressed for optimal use of the precious 
observing time at 10-m class telescopes). However, such biases 
do not seem capable of explaining two orders of magnitude of 
difference in the observed frequency of brown dwarf compan- 
ions with respect to planets and stellar companions. 

The dearth of close brown dwarf companions to solar-type 
stars, commonly referred to as the 'brown dwarf desert' (e.g, 
ICampbelletal .111 9881; iMarcv & Butleril2000l). i s seen to extend 
at wider separations. Only recently iPatel et al.l d2007l) have re- 
ported the first radial velocity detection of a handful of com- 
panions with M c sin i firmly established in the brown dwarf 
mass range at orbital distances of ~ 4 to 18 AU. These first 
results do not allow yet to produce an actual number for the 
brown dwarf frequency in this separation range. At larger sep- 
arations (~ 50 - 1000 AU), near-IR direct imaging surveys 
have confirmed a deficit in sub-stellar companions relative to 
stellar companions, but not quite as extreme as that apparent 
at orb ital radii within a few AUs. McCarthy & Zuckerman 
(2004) find / BD = 1 + 1 %, roughly a factor of ten lower 
than the stellar compa nion frequency in the same separation 
range. More recently Metc hev & Hillenbrand! (120091) have de- 
ri ved /bp = 3.2 + ^'%, a number f ormally compatible with that 
of McCarthy & Zuckerman (2004). In any case, the abundance 
of wide-separation brown dwarf companions is comparable to 
that of free-floating brown dwarfs relative to stars. The evidence 
for a not-so-dry desert at wide separations is usually interpreted 
as supporting the view that brown dwarfs form by core fragmen- 
tation just like stars. However, the finer details of the formation 
mechanism are not well understood, as mentioned above, and 
a more accurate characterization of the multiplicity of brown 
dwarfs is needed in order to shed light on the relative role of 
the various proposed scenarios. 

The dynamical determination of the mass of a few close 
brown dwarfs companions to low-mass stars and of brown dwarf 
binaries has been obtained for e clipsing systems for which 
both spectra can be observed (e.g., Zapatero Os orio et al.ll2004t 



Stass un et al.l 12006b . and by a combi nation of apertur e mask- 
ing interferometry and astro me try (llreland et al.ll2008h . As for 
the sample of companions to nearby dwarfs discovered by 
Doppler planet surveys with M c sin i values in the range be- 
tween high-mass planets and brown dwarfs, it is conceivable 
th at some of them a re in reality stars seen pole-on. Similarly 
to Halbwachs et al . (2000), several authors have attempted to 
combine the radial-velocity orbits with Hipparcos astrometry in 
order to determine the inclination a nd true mass o f the com- 
panions (for a review see for example ISozze tti 2009). Recently, 
Reffert & Ouirrenbach (2006) presented low-significance detec- 
tions of the astrometric orbits of the two outer companions 
in the HD 38529 and HD 168443 planetary systems, inferring 
masses in the brown dwarf regime of 37 My and 34 My for HD 
38529c and HD 168443c, r espectively. High-p recision astrom- 
etry with HST/FGS allowed Bean et all (120071) to determine an 
actual mass of 0.14 M e for the companion to HD 33636, origi- 
nally published with M c sin / = 9.3 My. Any successful attempt 
to resolve the ambiguity in the actual companion mass for such 
objects is thus worthwhile as it allows one to improve on the 
characterization of the multiplicity of brown dwarfs. This in turn 
permits to better understand the global nature of the sub-stellar 
companions found at larger separations by direct imaging sur- 
veys. 

In this paper we present a new attempt at combining the 
information from Doppler measurements with Hipparcos as- 
trometry to better the constraints on the mass of two Doppler- 
detec ted low-mass co mpanions to HP 1 31664 dMoutou et al.l 
120091) and HD 43848 dMinniti et al.ll2009l) . with published min- 
imum masses in the low-mass brown dwarf regime. A summary 
of the available data in the literature for the two systems is pre- 
sented in § [2] We derive improved contraints on the actual com- 
panion masses in § [3] Finally, in § |4]we a) put the new results 
in the context of the properties of the (few) close sub-stellar and 
massive planetary companions to nearby solar-type stars, col- 
lected in a catalogue table for ease of consultation and reference 
in future works, and b) discuss some of the implications of the 
present-day observational evidence for formation and structural 
models of massive planets and brown dwarfs. 

2. Doppler data and Hipparcos astrometry 

2.1. HD 131664 

Included in the Doppler search for southern ext rasolar planets 
carried out with the HARPS dPepe et all I2003I) spectrograph 
on the ESO 3.6-m telescope at La Silla Observatory, the bright 
G- type star HP 13166 4 (HIP 73408) was recently announced 
bv lMoutou et al.l (|2009) to be orbited by a companion with min- 
imum mass of M c sin/ = 18 My on an eccentric orbit of pe- 
riod P =s 5.3 yr (for convenience, see Table Q] for a summary 
of the properties of the primary and of the orbital parameters 
of detected companion). At the distance of HP 131664, the in- 
ferred orbital separation ranges between 35 mas and 100 mas (as 
pointed out by Moutou et al.), making it a potentially interest- 
ing target for future direct imaging observations. The viability 
of this investigation depends in particular on the true mass of 
the companion and the actual contrast ratio. In the former case, 
high-precision astrometr y with Gaia in space and VLTI/PRIMA 
from the ground (e.g., ISozzettil l2009, and references therein) 
will allow to derive accurate values of the actual mass of the 
companion. However, it must be noted that the minimum astro- 
metric signature induced on the primary is just under 1 mas, as 
opposed to a median single-measurement error cr HIP » 3 mas. 
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Table 1. Stellar characteristics of HD 131664 and HD 43848 and spectroscopic orbital elements for the two Doppler-detected 
companions. Reported errors are from the discovery papers. 
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HD 131664 


G3V 


8.13 


18.04 


55.43 


1.10 


1951 


52060 


0.638 


149.7 


359.5 


18.15 


3.17 












±0.03 


±41 


±41 


±0.02 


±1.0 


±22.3 


±0.35 





HD 43848 G2V 8.65 26.99 37.05 0.89 2371 53227 0.690 229.0 544 25.0 3.40 

±840 ±65 ±0.12 ±9.0 ±200 



Hipparcos observations of this star could then help place use- 
ful mass constraints on the companion mass. Indeed, the Double 
and Multiple Star Annex of the Hipparcos Catalogue reports a G 
flag, indicating that a 7-parameter solution (allowing for accel- 
eration in the proper motion) was found to significantly improve 
the standard 5-parameter single-star fit. The reported accelera- 
tion solution for HD 131664 is based on 96 datapoints, two of 
which corresponding to observations retained only by the NDAC 
consortium. 



2.2. HD 43848 

Doppl er measurements take n with the MIKE echelle spectro- 
graph dBernsteinetal.1 120031) on the 6 . 5-m Magellan II (Clay) 
telescope have revealed dMinniti et al.l 120091) the nearby solar- 
type star HD 43848 (HIP 29804) to be orbited by a companion 
with M c sin i = 25 M; on a high-eccentricity orbit with P « 6.5 
yr (for completeness, Table [TJ also reports the summary of the 
properties of the primary and of the orbital parameters of de- 
tected companion). This system also appears to be of potential 
interest for future direct imaging observations, with a separation 
at apoastron of « 0.2 arcsec. Its minimum astrometric signature, 
at the distance of HD 43848, is as 2.5 mas, a value compara- 
ble to the typical precision of Hipparcos astrometry for this star 
(median error cr HIP ~ 3.3 mas). An astrometric solution with 
acceleration terms for this star is also present in the Hipparcos 
Catalogue, based on a total of 76 measurements (with one ab- 
scissa retained only by FAST, one only by NDAC, and one re- 
jected by NDAC in the solution). For both HD 43848 and HD 
131664 there appears to be indication of the presence of a long- 
period trend in the Hipparcos data, rendering worthwhile a fur- 
ther investigation of the available astrometry. 



3. Combined radial velocity+astrometry solution 

When searching for evidence of the presence of orbital sig- 
nal in Hipparcos data of a given star due to a spectroscop- 
ically discovered low-mass companion, a typical procedure 
is applied in which information from radial velocities is as- 
sumed known, and one resorts to probe the region of the pa- 
rameter space not covered by spectroscopy . Orbital fits to the 
Hipparcos IAD are then usually performed (Mazeh et al. 1999; 
Halbwach set al.ll2000t jHan et al.ll2001b IZucker & Mazehfe OOO; 
Reffert & Ouir renbachll2006l: IWittenmver et al.ll2009l) by keep- 
ing fixed four orbital elements (P, e, To, ai) to their spectro- 
scopically determined values and by solving for inclination an- 



2 Note that Minn iti et all d2009h looked for, and failed to find, hints 
of higher dispersion in the Hipparcos measurements for HD 43848 with 
respect to those of stars of similar magnitude and distance. 



gle i and position angle of the ascending node Q, with the ad- 
ditional con straint that the astrometric semi-major axis satisfies 
the equality dPourbaix & J orissen 2000): 



a sin i = 9.19 X I0~ s PK Vl - e 2 7T+ mas, 



(1) 



where P is in days, the semi-amplitude of the radial velocity 
curve K is in m s , and the orbital parallax 7r+ is in mas. The 
resulting fitting procedure has then a total of 7 adjustable param- 
eters, i.e. the five astrometric parameters (two position offsets in 
right ascension and declination, two differential corrections to 
the proper motion in right ascension and declination, and the 
parallax) plus i and O. 

The experience of various authors indicates that one must 
be careful in not asking Hipparcos data to tell mo r e than they 
actu ally can. IPourbaix! (1200 II) . iPourbaix & Arenoul (|2001|) . and 
later IZucker & Mazehl (|2001) have for example shown that the 
Hipparcos IAD, while useful to put upper limits on the size of 
the astrometric perturbations, must be interpreted with great cau- 
tion when attempting to derive actual astrometric orbits for sub- 
stellar companions with semi-major axes close to or even be- 
low the typical single-measurement precision of the satellite. So 
while the astrometric orbit does not need to be actually detected 
in order to derive constraint s on the values of i and M c , when 
not outright refuted (e.g., IPourbaix & Arenoul 12001). reported 
detections usu ally have relatively low leve ls of statistical confi- 
dence (e.g., Reff ert & Ouirre nbach 2006). Based on the above 
considerations, we have adopted a multi-step approach to the 
treatment of the Hipparcos IAD for HD 131664 and HD 43848 
in an attempt to make a statistically solid case for any conclu- 
sions we derive. 

As a first step in the process, we decorrelated and weighted 
all available Hipparcos along-scan measurements for the two 
star s, following the prescriptions o f IPourbaix & JorisserJ ([2000) 
and Ivan Leeuwen & Evansl (119981) . We then asked whether an 
orbital model could result in an improved description of the 
Hipparcos data, and used the Th iele-Innes representation of 
a photocentric orbit dHeintz|["l978l) to carry out a linear least 
squares fit over a large grid of periods bracketing the ones ob- 
tained from the radial velocity measurements, while keeping 
fixed e and To to their spectroscopic values. We thus seeked to 
minimize 



x 2 = 



Av- 



(2) 



where the superscript t indicates transposed, Av is the vector of 
along-scan Hipparcos IAD residuals to a single-star, 5-parameter 
fit, dv/dpk is the vector of partial derivatives of the along-scan 
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Table 2. Orbital fits to the Hipparcos IAD for HD 131664 and HD 43848, assuming knowledge of P, e, T , and u from the 
spectroscopic orbit. Adjustable parameters were corrections to the five astrometric parameters, i, and £1 The resulting estimated 
masses of the companions are also given, based on the fitted value of the inclination angle. 



Star Name Aa AS h~fl~* AfT s i f2 M~c X$$) ~xW) P{F) 

(mas) (mas) (mas) (mas yr~') (masyr -1 ) (deg) (deg) (Mj) 



HD 131664 -2.8 ±0.6 -2.1 ±0.7 0.5 ± 0.9 -5.1 ±0.7 -4.1 ±0.8 55 ± 33 22 ± 28 23+f 1.39 1.15 0.001 



HD 43848 -4.9 ± 0.7 -9.8 ± 0.8 0.1 ±0.9 -8.6 ±0.3 0.3 ± 0.9 12 ±7 288 ± 22 120*'f 1.67 1.22 0.0002 



HD 131664 



CM 

X 




10 100 1000 10000 

Trial Period (d) 

HD 43848 




10 100 1000 10000 

Trial Period (d) 

Fig. 1. Periodograms of HD 131664 (top) and HD 43848 (bot- 
tom) assuming e and Tq from the known spectroscopic orbit. 



coordinate with respect to the k-th fitted parameter, Apt the rel- 
ative correction, and E -1 is the inverse of the covariance matrix 
of the solution. The fitted model is fully linear in 9 parameters, 
i.e. the five astrometric ones and the four Thiele-Innes constants 
A, B, F, and G. The results of this period search are shown for 
both stars in FigureQ] The presence of a long-period trend in the 
data appears clear for both HD 131664 and HD 43848, in clear 
accord with the choice of an acceleration solution for both ob- 
jects at the time of publication of the Hipparcos catalogue. The 
orbital peri od is well-known fro m spectroscopy in the case of 
HD 1 3 1 664 iMoutou et aLl (12009) report uncertainties on the or- 
de r of ~ 2%), but less so for HD 43848 (uncertainties reported 
bv lMinniti et all d2009l) are on the order of ~ 17%). Despite this, 



and the fact that Hipparcos only covered a fraction of the or- 
bit path of both objects, it is nonetheless worthwhile to investi- 
gate whether a full-fledged orbital solution can be extracted from 
the Hipparcos data, which robustly passes statistical screening. 
Following a procedure similar to that adopted by others and out- 
lined above, we then fitted the Hipparcos IAD for the two objects 
keeping P, e, and To fixed at their spectroscopic values, while ad- 
justing the five correction terms to the astrometric parameters 0, 
i and Q (with the additional constraint of Eq. Q}. We adopted a 
dense two-dimensional grid of starting values for both i and O, in 
order to solve a local non-linear minimization problem in which 
we used a model function expressed in terms of the Campbell 
orbital elements. In order to assess the impact of the more or 
less precise knowledge of the spectroscopic orbit, the procedure 
was re-run 100 times, each with a different set of (fixed) spectro- 
scopic orbital elements drawn from gaussian distributions with 
mean and standard deviation the best-fit value and its formal er- 
ror, respectively. 

The results of the orbital fit to the Hipparcos IAD of both HD 
131664 and HD 43848 are summarized in Table [2] The reported 
uncertainties on i and Q, as well as the derived mass values for 
the companions, take into account the 1-cr formal errors on the 
parameters of the spectroscopic orbits. The inferred true masses 
for HD 131664b and HD 43848b are 23+f Mj and 120^f My, 
respectively. Taken at face value, these numbers would imply 
that the companion to HD 131 664 has a mass in the brown dwarf 
regime at the 2 - cr confidence level, while the unseen object 
around HD 43848 is likely a low-mass M dwarf (at the 1-cr 
confidence level). An attempt at using the FAST and NDAC data 
separately for both stars resulted in solutions formally compati- 
ble with the ones reported here, albeit with larger uncertainties 
and looser constraints on the derived mass estimates. 

In order to assess the statistical soundness of the derived as- 
trometric orbits, several checks can be performed. First, based 
on an F-test of the null hypothesis that there is no companion, 
we find that the addition of two parameters to the model de- 
scribing the Hipparcos IAD significantly improves the fit in both 
cases: P(F) = 0.0002 and P(F) = 0.001 forHD 131664 and HP 
43848, respectively. Second, similarly to Reffert & Ouirrenbach 
d2006l) . we have further explored the reliability of the orbital so- 
lutions by inspecting the joint confidence regions (corresponding 
to iso-^ 2 contours) in the i - Q plane, for a given statistical level 
of confidence. As shown in Figure|2] the 95% (2-cr) iso-^- 2 con- 
tours cover relatively narrow regions in the i - Q plane, ruling 
out with confidence very small inclinations in the case of HD 



3 We used in the fitting proc ess both and its transf ormed value 
n '* = log 71 ** a trick utilized by Pourbaix & Jorissen (2000) to prevent 
parallaxes from becoming negative, and converting back to 7r* at the 
end of the minimization procedure. No significant differences in the 
final results were found. 
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Fig. 2. Iso-^ 2 contours corresponding to confidence intervals in 
two dimensions containing 68.3% (solid line), 90% (dashed- 
dotted line), and 95% (dashed line) of the values of i and Q. as 
obtained in a 7-parameter fit to the Hipparcos IAD of HD 43848 
(top) and HD 131664 (bottom), assuming P, e, To, and to from 
the known spectroscopic orbit. The best-fit (i,Cl) solutions for 
both stars are indicated with stars. 
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Fig. 3. Distribution of the ratio a/cr a of false astrometric or- 
bits for single stars observed by Hipparcos, approxim ated by a 
Rayleigh-Rice law following Hal bwachs et al. (2000). The two 
arrows indicate the values of a/cr a obtained for HD 43848 and 
HD 131664 based on a Monte Carlo bootstrap simulation. 



131664, and close to edge-on configurations in the case of HD 
43848. Ambiguities in the orbit orientation are clearly seen in 
both cases, highlighted by local minima in the x 1 surface corre- 
sponding to configurations with the opposite sense of revolution. 
Nevertheless, one would conclude that Hipparcos astrometry has 
successfully detected orbital motion induced by the two mas- 
sive companions to HD 131664 and HD 43848. We consider the 
95% confidence level as the appropriate measure of the degree of 
statistical robustness with which these results can be reg arded. 
Final ly, we have applied a standard bootstrap method (Efron 
119821) . based on 1000 Monte Carlo resamples with replacement, 
to the available data for both stars in order to derive empirical er- 
ror estimates on the derived orbital semi-major axes (constrained 
by Eq.Q] du ring the fits). Our proced ure closely follows the one 
adopted by Zucker & Mazeh (2001), who first realized, due to 
the assigned correlation between pairs of Hipparcos measure- 
ments for which both FAST and NDAC data are available, the 
need to preserve the pairing of the measurements while drawing 
new random datasets, in order to make the bootstrap method ap- 
plicable in the first place. The two resulting ratios a/cr a have 
then been compared to the Rayleigh-Rice law distribution of 
false a/cr a one can exp ect to obtain in the cas e of Hipparcos 
data of single stars (see (Halbwachs et al. 2000)). As shown in 
Figure [3] the nominal values of ajcr a for HD 131664 and HD 
43848 fall in the tail of the Rayleigh-Rice distribution: the prob- 
ability that a/<r a exceeds the observed one is found to be 0.03 
and 2 x 10~ 5 for HD 43848 and HD 131664, respectively. This 
can be again interpreted as a detection of a significant Hipparcos 
astrometric orbit for both stars, at the 97% and 99.99% confi- 
dence level, respectively. For the purpose of this study, we adopt 
the more conservative above mentioned 95% confidence level to 
gauge the actual degree of statistical confidence with which the 
orbits of HD 43848b and HD13664b are considered detected. 

4. Summary, discussion, and conclusions 

We have inspected the Hipparcos IAD for two stars, HD 43848 
and HD 131664, with Doppler-detected companions with min- 
imum masses in the brown dwarf regime. We have presented a 
body of supporting evidence which appears to confirm the de- 
tection of orbital motion in the Hipparcos IAD of HD 43848 and 
HD 131664, at a somewhat significant (95%) level of statistical 
confidence. The inferred actual masses of HD 43848b and HD 
131664b are found to be M c = 120^ 7 Mj and M c = 23+f 
Mj, respectively. The former thus appears to be a late M dwarf 
(at the ~ 1 — a confidence level), while the latter appears to 
be a brown dwarf (at the ~ 2 - cr confidence level). Taken at 
face value, HD 131664b is nominally the lowest-mass brown 
dwarf confirmed with a combination of Hipparcos and precision 
Doppler measurements (the astrometric orbits of giant planets 
obtained with HST/FGS + RV data lie in a different ballpark). 
The larger uncertainties in the spectroscopic orbital elements of 
HD 43848, combined with the low value of i inferred from the 
orbital fit, translate in larger uncertainties in the mass estimate, 
leaving margin for this unseen object to also be a brown dwarf 
(as well as a higher-mass M dwarf). 

Ba sed on the mass-luminosity relations of Del fosse et al.l 
(120001) , HD 43848b, with its nominal mass of 120 Mj, is ex- 
pected to have magnitude differences of about 8.7, 5.7, 5.4, 
and 5.2 in the V, J, H, and K bands, respectively. With a pro- 
jected separation of about 0.15 arcsec at apoastron it should 
be detectable with current AO instruments at 8m class tele- 
scopes. The magnitude difference between HD 131664b and 
its parent star depends critically on the age of the system. 
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iMoutou et ail d2009h give an age of 2.4 + 1.8 Gyr for HD 
131664 based on isochrone fitting. Additional clues can b e de- 
rived from stellar activity indicators. IMoutou et al.l (12009) mea- 
sured logR HK = -4.82 ± 0.07. Thi s corresponds to an age 
of 3.5 Gyr using the calibration by iMamaiek & Hillenbrand! 
(2008). We also searched for X-ra y emission from t he system. 
The ROSAT Faint Source Catalog dVoges et al.ll2000l) includes a 
source (1RXS J150003.3-733128) at 40 arcsecfromHD 131664 
(with a quoted positional error of 29 arcsec). The association is 
then doubtful. Assuming 1RXS J150003.3-733128 is the X-ray 
counte rpart of HD 131664, we derived (using the flux calibra- 
tion bv lfflinsch etal.1 (1 19991) ) log L x = 28.5 and an age of 2 .1 
Gyr (using the calibration by Mamaie k & Hillenbrand! ([2008)). 
Alternatively, a null detection would imply an older age. Stellar 
activity indicators then exclude the youngest stellar ages com- 
patible with isochrone fitting. A more plausible lower limit to 
the stellar age is about 1.5 Gyr. HD 131664b, given our best-fit 
mass value and the age of 2.4 Gyr, is expected to have magnitude 
differences of about 15.0, 15.2, and 17.2 in /, H, a nd K bands 
respectively, based on the models of Bara ffe et al.l (12003 ). The 
corresponding effective temperature is about 700 K, at the cool 
end of the currently known T dwarfs. Considering the small pro- 
jected separation (< 0.1 arcsec), such a contrast is not achievable 
with current instrumentation, and it would also be challenging 
for the next generation of direct imaging instruments. 

Any successful attempt, such as the one presented here, 
at deriving true masses of the substellar companion candi- 
dates detected from radial velocity surveys is definitely worth- 
while, as the results can be seen in the context of the ob- 
served paucity of close brown dwarf companions to solar-type 
stars in the solar neighborhood, the well-known brown dwarf 
desert, and in particular their binarity properties. For example, 
taken at face value the rather large eccentricity of HD 131664b 
(e = 0.638) supports the notion that massive planetary com- 
panion s and brown dwarfs are prefer ably found on eccentric 
orbits dRibas & Miralda-Escu de 2007). On the other hand, the 
high metallicity of the parent star ([Fe/H=+0.32) does not follow 
the prop osed trend of metallicity vs. m ass for substellar com- 
panions (Ribas & Miralda-Escude 2007). While not conclusive, 
these findings bring new/updated information that can certainly 
help to better our understanding of the formation mechanism of 
such objects. 

4.1. Properties of massive planets and brown dwarf 
companions to solar-type stars 

To put our discussion in a more general context, we summarize 
in Table [3]the main properties of the companions with projected 
masses between 10 and 80 My, and semimajor axis smaller than 
~ 7 AU, orbiting main-sequence stars with masses between 0.7 
to 1.5 Mq. This selection matches the sample of most high- 
precision radial velocity surveys and avoids the additional com- 
plications of large variations in the stellar mass and evolution- 
ary status (very small number statistics being one of the most 
relevant). The literature data collected here are meant to pro- 
vide ease of consultation and reference for future works on the 
many outstanding issues we touch upon thereafter. In the Table, 
the determination of or constraints on the inclination and true 
companion mass derived from astrometry or transit photometry 
are included when available. For completeness and reference, 
we also include companions with substellar projected mass and 
astrometrically derived stellar mass. We have also included HD 
33636, whose companion with M c sin i just below the adopte d 
threshold was shown to be a low mass star dBean et al.l 12007). 
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Fig. 4. Period vs eccentricity (top panel) and eccentricity vs 
companion mass (bottom panel). In the lower panel plot only 
objects with period longer than 20 days are shown). Red filled 
circles: companions with true estimated masses in the range 10- 
80 Mj. Blue triangles: companions with projected masses in 
the range 10-80 Mj. Green squares: companions with projected 
masses in the range < 80My but true stellar masses > 80M/. The 
three objects with larger symbols belong to the multiple systems 
HD 38529, HD 168443 and HD 202206. 



Orbital elements and pro jected m asses are from l i sted d iscovery 
papers, iButler et al.l d2006l) or lHalbwachs et alj (l2000h. S tellar 
masses and metallicities are from lValenti & Fischerl d2005). dis- 
covery papers or additional references. Further information, such 
as additional stellar or planetary companions, the sources for the 
metallicity values, and the specific treatment of some of the mass 

and error estimates, can be found in th e notes. 

As discussed in the Introduction, Halbwachs et al. ( 2000) 
have shown how a significant fraction of Doppler-detected can- 
didate substellar companions were in fact low-mass stars viewed 
at low inclination. Nevertheless, a few candidates have masses 
firmly in the substellar regime, partially filling the brown dwarf 
desert. Interestingly, in two cases (HD 38529 and HD 168443) 
additional companions in the planetary mass regime were found. 
Q We note that such system configurations are more typical of 
planetary systems than of multiple stellar systems. 



4 Two additional probable cases of a system with a brown dwar f and 
a planet in c lose orbits are those of HAT-P-13 (Bakos et al. 2003) ar, d 
HD 202206 dCorreia et alj|2005l) . In the first case, the planetary nature 
of the inner companion is confirmed by the occurrence of the transit. 
The true masses of the more massive companions are not available but 
there are perspectives for deriving them thanks to the strong dynamical 
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Fig. 5. Orbital period vs companion mass (top), metallicity vs 
companion mass (center), stellar mass vs companion mass (bot- 
tom). Symbols as in the previous Figure. 

In Figure H and Figure [5] we show a set of correlation di- 
agrams among the most relevant quantities of the systems col- 
lec ted in Tabled 

iRib as & Miralda-Escudl d2007l) noted a possible trend in 
the eccentricity-mass distribution of giant planets, with more 
massive planets (M c sin/ > 4 My) orbiting on average with 
larger eccentricities. The recent analysis by Wright et al. 
(200 9), who used 1 My as the cut-off value, confirm this trend. 
The IRibas & Miralda-Escudel d2007l) analysis also suggested 
that the metallicity of stars hosting massive planets and brown 




0.4 0.6 

Eccentricity 




interactions in the case of HD 202206 and the transit timing variations 
in the case of HAT-P-13. 



[Fe/H] 

Fig. 6. Upper panel: cumulative distribution of the eccentric- 
ity of objects with period longer than 20 days orbiting main- 
sequence stars with mass between 0.7 and 1.5 M H . Continuous 
(black ) line: planets with M c sinj < 4 My (from iButler et al.1 
(l2006h). Dotted (purple) line: planets with 4 < M c sin < 10 My 
(from But ler et all (|2006)). Dashed (red) line: objects with mass 
10 < M c < 80 My (this paper), Dashed-dot (blue) line: objects 
with 10 < M c sin; < 80 My, but excluding objects with astro- 
metric masses in the stellar regime (this paper). Lower panel: 
cumulative [Fe/H] (with no period cut). Same symbols as in the 
upper panel. 



dwarf candidates is c ompatible to that of the solar neighborhood 
([Fe/H]^ -0.15, e.g. lNordstrom et al.ll2004l) . differing at the 3-cr 
level from the super-solar metallicity distribu tion of stars host- 
ing lo wer mass planets ([Fe/H]^ +0.24, e.g. iFischer & Valentil 
2005). Hints of su ch a trend had alr eady been seen in a more 
limited dataset by iRice et aL (2003). We revisit here this sta- 
tistical analysis, considering all objects included in Table [3] as 
well as t he sample of exop lanets with M c sin/ < 10 My re- 
ported in Butl er et ail d2006l) . For this purpose, we have divided 
the full sample in three subsets according to minimum mass: 
M c sin /< 4 My, 4 < M c sin i < 10 My, and 10 < M c sin i < 80 
My, plus a fourth subsample formed by those objects in Table [3] 
with estimated true masses in the range 10 < M e < 80 My. We 
have then performed on these subsets a Kolmogorov-Smirnov 
(K-S) test, to measure to what extent their e and [Fe/H] distri- 
butions might differ, and a Wilcoxon Rank-Sum (R-S) test, to 
measure to what degree the mean e and [Fe/H] might be con- 
sidered different. The resulting values of the probabilities (Pr(D) 
and Pr(Z), respectively) of the null hypothesis for both tests on 
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Table 4. Results of the K-S and R-S tests on different subsets 
of systems: all planets from the Butler et al. (2006) catalog with 
M c sin 2' < 4 Mj (1), companions with 4 < M c sin i < 10 My (2), 
companions with 10 < M c sin/ < 80 My but excluding those 
with confirmed mass larger than 80 My (3), and companions with 
true masses in the range 10 < M c < 80 My (4) . 



Sub-sample 


D 


Pr(D) 


Z 


Pr(Z) 


Eccentricity Distribution 


1 vs 2 


0.435 


0.002 


3.142 


0.001 


1 vs 4 


0.50 


0.08 


2.18 


0.01 


1 vs 3+4 


0.449 


0.001 


3.4535 


0.0003 


2 vs 4 


0.32 


0.64 


0.73 


0.23 


2 vs 3+4 


0.21 


0.65 


0.24 


0.41 


Metallicity Distribution 


1 vs 2 


0.22 


0.27 


-0.67 


0.25 


1 vs 4 


0.46 


0.06 


-1.57 


0.06 


1 vs 3+4 


0.389 


0.002 


-2.745 


0.007 


2 vs 4 


0.31 


0.51 


-1.11 


0.13 


2 vs 3+4 


0.28 


0.23 


-1.36 


0.09 



the various subsamples are reported in Table H] while Figure [6] 
shows a comparison between the resulting cumulative distribu- 
tions of e and [Fe/H] for the abovementioned subsets of systems. 

From investigation of Table [4] and Figure [6] a few conclu- 
sions can be drawn. As for the e distribution, the results further 
corroborate the notion that eccentricities of not so massive plan- 
ets are clearly less pronounced than those of more massive plan- 
ets and brown dwarf companions^. On the other hand, massive 
planets appear to have an overall e distribution which is indis- 
tinguishable from that of brown dwarf candidates. Similarly, the 
[Fe/H] distribution of the hosts appears significantly different for 
companions with M c sin /< 4 My and for brown dwarf compan- 
ions, while other trends are marginal, and in particular metallic - 
ities of massive planet hosts and of primaries with brown dwarf 
companions are statistically the sameQ 



4.2. Implications for formation and structural models of 
massive planets and brown dwarfs 

The above evidence can be interpreted, keeping in mind the rel- 
atively small-number statistics regime we're dealing with here, 
in the context of the proposed formation scenarios and internal 
structure models of high-mass planets and brown dwarfs, and as 
a function of the main properties of the stellar hosts (binarity, 
mass, metallicity). 



5 The statistical analysis presented here does not take into account 
the eccentricity bias in Keplerian fits to Doppler data collected by 
radial-velocity planet search programs, which a) underestimates the 
abundance of low-amplitude, low-eccentricity planets, and b) makes 
low-amplitude, high-eccentricity planets hard to unveil ( Shen & Turner 
2008; lQ'Toole et alJl2009l : IValenti et al.ll2009l) 

6 It also appears that the brown dwarfs with the largest actual mass 
estimates seem to be found more frequently as close companions to 
stars belonging to the thick disk population of the Milky Way (see 
TableO, but this can only be regarded as a tentative speculation, given 
the very limited statistics, large mass uncertainties, and heterogeneity 
of the datasets. 



4.2.1 . Formation and orbital evolution of massive planets and 
brown dwarfs 

As already remarked in the Introduction, the lower mass limit 
for the formation of self-gravitating objects from fragmenta- 
tion of molecular cloud cores is today though t to be on the 
order o f a few Mi (e.g.. IWhitworth etal1l2007l: ILuhman et ail 
20071). iRibas & Miralda-Escudel d2007l) and, more recently, 
Font-Ribera et alj (120091) proposed direct cloud fragmenta- 
tion, followed by inward migration by disk capture, as a 
formation mechanism for the high-mass tail of the plane- 
tary population. The detection of a significant number of 
free-floating objects down to ~ 6 My and the observed 
continuity of the subs tellar mass function in young star 
forming regions (e.g.. ICaballero et al.l 120071 and references 
therein; Zuckerman & Songl 120091 and references therein) are 
indications that indeed a star-like formation process may form 
brown dwarfs as w ell as planetary mass ob jects Q. Indeed, the 
process outlined by Font-Rib era et al.l {2009} should be effective 
independently of the metallicity of the parent cloud, and is ex- 
pected to produce an eccentricity distribution similar to that of 
binary systems. The results shown in Tableland in Figure[6]can 
be seen as supportive of this scenar io. On the other hand, th e disk 
capture mechanism proposed by iFont-Ribera et all J2009) does 
not naturally explain the existence of the brown dwarf desert, as 
in their model increasingly more massive objects should actu- 
ally be easier to capture. Also, the observed differences in the 
mass distribution between isolated obje cts and companions to 
solar-t ype stars, the latter, as shown by Greth er~& Lineweaverl 
(120061) . exhibiting a minimum (the "driest part of the brown 
dwarf desert") at « 30 My (31^ My), do not appear to be read- 
ily explainable by this mechanism. 

Massive planets and brown dwarfs can also be formed 
by gravitational fragmentation of ext ended disks around solar- 
type p rimaries. In the models of e.g. Stam atellos & Whitworfhl 
(2009), brown dwarfs form by direct gravitational collapse only 
in the outer regions of the disks, and tend to be scattered further 
out, or even into the field. This mechanism can explain the ex- 
istence of the brown dwarf desert at small separations. The two 
other conclusions that can be drawn from these models are that 
1) no planets formed by disk instability can really be found as 
close companions, as they can only form in the outer regions 
of the disks and are scattered out with much higher efficiency 
than brown dwarfs, and b) as a direct consequence short-period 
Doppler-detected massive planets and brown dwarfs should be 
preferentially low mass stars. On the one hand, this mechanism 
can be seen as supported by the data presented here (the proper- 
ties of massive planets, brown dwarfs and low mass stellar com- 
panions being very similar). On the other hand, some of the ob- 
jects in the sample of Table [3] can hardly be explained by this 
approach (e.g., the transiting objects XO-3b and Co RoT-3b). 

In the sti l l the o retically debated (e.g., iMaver et al. l 

l200l iRafikovl 12001 IStamatellos & Whitwortfl 120081: iBossi 
120091 for a review see iDurisen et al.l 120071) disk instability 
model of giant planet formation, massive planets with relatively 
eccentric orbits and moder ate to large orbita l separations a re its 
more likely product (e.g., iRice et all 120031: lBolevll2009l) . No 
clear prediction of the expected shapes of the orbital elements 
and mass distributions of planets formed by disk instability 
is yet available, due to the numerical complexities of the 
simulations and to some still open theoretical issues on the input 



7 However, see e.g. IBossi fcOOd) and iBate et al.l d2002h for alterna- 
tive scenarios in which isolated planetary mass objects are the result of 
ejections through dynamical interactions in protoplanetary disks 
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physics (see e.g. Boss 2009, and references therein). However, 
if high-precision measurements of the actual masses of massive 
planets as well as of those of the substellar companions listed in 
Table [3] were to become available, and they were to prove that 
these objects are not preferentially low mass stars, this evidence 
would argue in favor of a common formation mode for massive 
planets and brown dwarfs. 

In the more widely accepted core-accre tion mecha- 
nism f or the forma t ion o f giant planets (e. g., iPollack et alj 
[19961 lAiibert et alj l200l llda&Linl l200l for a review 
see iLissauer & Stevensonll2007h very massive planets and low- 
mass brown dwarfs on e ccentric orbits are not a natural out- 
come (e.g., IPollack et al.1 



120051: iKennedv & Kenvon 



1996t [Alibert et al. 2005; llda & Linl 



2008). However, recent models have 



shown that this mechanism might also be capa ble of forming 
such m assive objects. Indeed, in the model of iMordasini et al] 
(2009) planets with masses as large as 38 Mj are formed 
in long-lived massive disks around a solar mass star. Given 
the increasingly lower probability of forming very massive 
planets or brown dwarfs ('d euterium-burning' planets in the 
words of Bara ffe et al.l d2008|)). the brown d warf desert appears 
compatible with the Mordasini et al. (2009) model (for exam- 
ple, the probability of forming objects with masses exceed- 
ing the deuterium- burning thresho l d is f ound to be ~ 0.4%, 
in accord with thelMarcv & Butler! d2000t) estimate). While the 
Mordasin i et alj d2009t) model does not include multi -planet sys- 
tems and planet-planet interactions, likely to play an important 
role in shaping the observed mass, period and eccentricity distri- 
butions, we note that the expected location of such very massive 
planets is between 2 to 5 AU, similar to that observed for HD 
38529c, HD 168443c and the object studied in this paper, HD 
131664b. However, no prediction of this model on the eccentric- 
ity distribution is available, which would be a critical element for 
discriminating between this mechanism and the other aforemen- 
tioned options. In addition, super-planets and low-mass brown 
dwarfs do exist also at sm all separation, where the y are not ex- 
pected on the basis of the IMordasini et ail (2009) model (e.g., 
the t ransiting systems HAT-P-2b dBakos et al.l 120071) . WASP- 
14b dJoshi et al.[|2009h . XO-3b, and CoRoT-3b, but a lso other 
massive o bjects such as HP 1 62020b dUdrv et al.ll2002l) and HD 
41004Bb dZucker et alj|2004 ). 

Finally, of particular interest are the initial claims of a pos- 
sible correlation between massive planets, eccentric orbits, and 
large values of the angle between a planetary orbit and the stel- 
lar rotation axis, as determined from spectros copic measure- 
ments of the Rossi ter-McLaughlin (R-M) effect dRossiterill924l; 
iMcLaughlinl 1 1924b in transiting systems. The XO-3, WASP- 
14, and HD 80606 systems all have a close-in, massive planet 
on a very eccent ric orbit and with a s ignificant spin-orbit mis- 
alignment (e.g., iJohnson et al.l H009, and references therein) 
0. At first glance, this might imply that the orbital migration 
history of massive, eccentric exoplanets is somehow different 
from that of less massive close-in Jupiters, However, the pic- 
ture is likely not to be so simple, not only because of the evi- 
dence for other transiting systems with massive, eccentric plan- 
ets (HAT-P-2, HD 17156) showing no signs of spin-o rbit mis- 
alignment (e.g.. lWinn et al.|[2007t Barbieri et al. 2009), but also 
for the recent measurements of the R-M effect in transiting 



8 While the results are not conclusive due to the presence of systemat- 
ics in the dataset, it is worth noting tha t the CoRot-3 syste m also exhibits 
a formally non-zero spin-orbit angle (Triaud et al. 2009). However the 
large uncertainties reported make the detection of spin-orbit misalign- 
ment in the CoRoT-3 system only marginal. 



systems with Jupiter-mass planets o n circular orbits , such as 
CoRot-1, HAT-P-7, and WASP-17 ( iPont et al.l 120091 and ref- 
erences therein) which point at large values of the projected 
spin-orbit angles. Overall, the evidence collected so far is a 
likely indication of the variety of possible outcomes of the 
complex process of migration to close-in orbits of compan- 
ions with a wide range of masses, which include a host of 
proposed dynamical mechanisms (e.g., interactions between a 
planet and the gaseous/planetesimal disk, planet-planet resonant 
interactions, close encounters between planets, and secular in- 
teractions with a companion star) as well as different forma- 
tio n scenarios. Statistical s tudies such as the one carried out 
bv lFabrvckv & Winn] d2009l) . who showed the emergency of a bi- 
modal distribution of spin-orbit angles, on an increasingly larger 
sample of transiting systems will ultimately be the optimal way 
to compare an ensemble of measurements of the R-M effect with 
the predictions of migration theories. 



4.2.2. The impact of binarity 

Among the distinctive features of exoplanets discovered around 
mem bers of various types o f binary and triple systems 
(e.g., Eggenberger & Udry 2009, and references therein), one of 
the most intriguing is the evidence that solar-type stars members 
of multiple stellar systems appear to be pre ferential hosts of the 
most massive planets in short-period orbits (Desidera & B arbieril 
2007), and the fact that the planetary companions with the high- 
est eccentricities all h ave a stellar or brown dwarf compan- 
ion dTamuz et al.l2 008). These trends seem to indicate that planet 
formation and/or migration in binaries may proceed differ- 
ently than around single stars. Indeed, theoretical studie s (iKleyl 
2000; IWu & Murravll2003l; iFabrvckv & Tremaindl2007l) within 
the context of the core accretion model of giant planet for- 
mation suggest that the presence of a fairly close companion 
significantly enhances the growth rate and make the migration 
timescale of the planet shorter. On the other hand, the massive 
super-planets or low-mass brown dwarfs found at separations of 
about 2-3 AU orbit stars that are single or with very wide com- 
panions (HD 38529 has a 0.5 M Q common proper motion com- 
panion at a projected separation of 12000 AU), compatible with 
long-lived , undisturbed disks req uired to form such objects ac- 
cording to lMordasini et al.ld2009l) . Models of giant planet forma- 
tion by disk instability, however, come to opposite conclusions, 
with giant planet formation s ignificantly suppre ssed in binaries 
with separations < 100 AU (Ma yer et al.ll2005h . From inspec- 
tion of Table [3] the abovementioned trends seem to be supported 
only in part. For example, the short-period objects HD 98230b 
and HD 283750b orbit one of the components of wide bina- 
ries, but this does not appear to be the case for HD 162020b. 
Investigations are encouraged which would aim at verifying the 
possible existence of binary companions to CoRoT-3b and XO- 
3b. Furthermore, the brown-dwarf candidates with the highest 
eccentricities do not seem to be preferentially found in multiple 
systems. 



4.2.3. The role of the primary mass and metallicity 

In the core-accretion model of giant planet formation, the up- 
per limit on the mass of the planetary companion and the final 
orbital arrangement are expected to depend on stellar mass and 
metallicity. Based on arguments of protoplanetary disk size and 
lifetime as a function of M+, one would expect massive planets 
to be found with higher probabilities around more massive pri- 



10 



A. Sozzetti & S. Desidera: Hipparcos preliminary astrometric masses for the companions to HD 131664 and HD 43848 



maries dKennedv & Kenvon 2008 ), and at typically moderate to 
large separations dBurkert & Idal2007tlCurriel2009l) . Indeed, the 
observational evidence indicates that, for larger stellar masses, 
massive companions with > 10My are si gnificantly more nu- 
merous than ar ound solar type st ars (e.gjLovis & Mavorj l2007; 
|Johnsonir2008i For a review see lHatzesI d2008l) . and references 
therein). On the opposit e end, M dwarfs are expected to show a 
paucity of giant planets (Laughlin et al. 20 041: llda & Linl | 2005l) . 
which is also obs erved ( Endl et al. 2006t iJohnson et alj 12007; 
Baile y et al.l2009l) . The core accretion mechanism also naturally 
predicts that super-massive planets should not be found at all 
around m etal-poor stars, g iven the lack of material for accumula- 
tion (e.g. . llda & Linl2005l) . The alternative disk instability mech- 
anism is instead rather in sensitive to t he val ues of M* and [Fe/H] 
of the stellar host (e.g., lBossll2002l [2006). As a consequence, 
a qualitative prediction of this model (e.g., iRice et aiTl2003l) 



is that massive planets found on moderate to large separations 
around massive and/or metal-deficient hosts would be likely to 
have been formed by gravitational instability. In this respect, 
the data collected here and the results of the statistical analysis , 
that corroborate the findings of Ribas & Miralda-Escude ( 2007), 
can be read as partly supportive of the latter view. Interestingly, 
the frequency of (massive) planets (f p ) around intermediate- 
mass ( mostly giant) stars appears to be rather independent of 
[Fe/H] dHatzesll2008l) . However, the four known planet - brown 
dwarf systems (HD 38529, HD 168443, HD 202206, HAT-P-13) 
all have super-solar mass, very metal-rich (main-sequence) pri- 
maries. The latter evidence might be seen as supportive of the 
core accretion formation mode, while the former dataset might 
point to a significant role of the disk instability mechanism. 
However, the global picture is likely to rather complex. For ex- 
ample, Doppler surveys for giant planets around intermediate- 
mass stars typically include targets which are evolved to some 
degree, given that massive main sequence stars are unsuitable 
for high-precision radial-velocity measurements (too few spec- 
tral lines, often broadened by high rotation rates). The very dif- 
ferent dependence of f p on [Fe/H] for intermediate-mass stars 
with respect to their solar-mass counter-parts might then reflect 
a non-primordial or igin of the metallic ity enhancement in solar- 
type planet hosts (Pasq uini et alj |2007) rather than point to dif- 
ferent formation modes. 

As for the possible dependence of brown dwarf frequency 
on stellar mass and metallicity, this is still a poorly understood 
issue. A quick look at literature data allows to speculate on 
the possibility that the brown dwarf desert may not be so 'dry' 
when it comes to close sub-stellar companions to intermediate- 
mass stars, in light of a handful of systems containing at least 
one c ompanion with min imum mass in the brown dwarf regime 
(e.g.. lOmiva et alj 12009, and references therein). In particular, 
the rec ently announced system of brown dwarfs around BD+20 
2457 dNiedzielski et al.ll2009l) has an architecture very similar 
to that of HD 168443, with a solar-type primary. It is reminis- 
cent of an origin in a massive circumstellar disk, further sug- 
gesting a scenario in which more massive sub-stellar compan- 
ions are found around more massive stellar hosts (and their fre- 
quency also increases with increasing primary mass). Such sys- 
tems would be more likely to have form ed by local gravit ational 
instabilities in protoplanetary disks (e.g jRice et al.ll2003l) rather 
than protostellar cloud fragmentation, given the difficulties in 
formi ng extreme mass-ratio binaries by the latter mechanism 
(e.g., Bate 2000). However, the combination of small-number 
statistics, the different priority given to observations of stars with 
very massive companions in Doppler surveys, the variable detec- 
tion thresholds as a function of companion mass and separation 



due to the decreasing radial-velocity precision in increasingly 
higher mass stars, and the uncertainty on the actual mass values 
for most of the objects included in Table[3](as well as those ref- 
erenced above) prevents one, as of today, from drawing any se- 
rious conclusions on the possibility that the brown dwarf desert 
may move in mass and/or separation range depending on M* and 
[Fe/H]. 

4.2.4. Structural and atmospheric models 

Finally, one of the most effective ways of distinguishing be- 
tween massive planets and brown dwarfs is through a compar- 
ison of their internal structure properties and of th e composi- 
tion o f their atmospheres (e.g.. IChabrier et aTl 120091 [BurgasseJ 
2009). Studies at this level are still in their infancy status, 
particularly for what concerns the possibility of determining 
the actual nature of such objects based on their condensate 
cloud formation prop erties, non-equili brium chemistry, and at- 
mospheric dynamics (Burgasser 2009). Howewer, the class of 
transiting systems is already pro viding data of relevance. For 
example, Mor dasini et alj Q2009) have proposed that massive 
planets might have a large content of heavy elements (about 
0.8 Mj for the individual case shown in their Fig. 12). In the 
case of the transiting massive planet HD 14 7506b (HAT-P-2b) , 
with a mass of 9.04 My and R = 0.982 Rj dBakos et al.ll2007D . 
just below the threshold of 10 Mj a dopted here fo r incl usion 
in Table [3] iBaraffe et all d2008l) and iLeconte et all d2009l) ar- 
gue from the measured radius for a total amount of heavy el- 
ements of about lMj, excluding a gaseous H/He object with 
solar composition. This indicates that indeed 10 Mj planets 
might be formed by the core accretion mechanism, unless such 
a l arge mass is the res ult of planetary collision, as speculated 
bv lBaraffe etaT] d2008l) . In the case of the super-massive transit- 
ing planet XO-3b, despite its rather uncertain radius (due to dis- 
crepancies in the stellar radius estimates from photometry and 
spectroscopy), there seems to be no need for a large content of 
heavy elements in the core dWinn et alj 120081)). Simila r ly, the 
22-My transiting object CoRoT-3b fits the IBaraffe et al J £2008) 
models fo r solar composition, without need for large metal en- 
richment dLeconte et alj|2009l) . Both HAT-P-2b and CoRoT-3b 
have surface gravities similar to those of young, very low-mas s 
brown dwarfs such as AB Pic (see Fig ure 2 of lBurgasserll 2009). 
suggesting an overlap in the parameter space of some physical 
properties (e.g., gas pressure, temperature). Interestingly, as dis- 
cussed by Deleuil et alj d2008l) . all transiting massive planets or- 
bit stars more massive than the Sun, as predicted by the core 
accretion model of planet formation. On the other hand, their 
actual existence on such short periods and relatively eccentric 
orbits (except for CoRoT-3b), and the fact that the metallicities 
of their parent stars are not skewed towards super-solar values 
calls for the possibility of differences in their origin. 

In conclusion, the present-day evidence on the orbital and 
structural properties of massive planets and brown dwarfs as 
close companions to nearby solar-type stars can be interpreted 
as supportive of a picture in which different mechanisms for 
the formation of such objects are at work. The above discus- 
sion also indicates how the lack of a clean, statistically signifi- 
cant sample of high-mass planets and brown dwarf companions 
to solar-type stars with well-determined mass estimates hampers 
the possibility to conclusively discriminate among the various 
competing modes of formation of such objects. From an ob- 
servational viewpoint, improvements in the determination of the 
multiplicity properties of sub-stellar companions are very much 
needed, and they will be obtained in the near future through 
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the combined contribution of a variety of techniques. For ex- 
ample, some of the massive planets and brown dwarfs objects 
in Table [3] are within reach of n ext generation dire ct imag- 
ing instruments such as SPHERE dBeuzit et al.l l2QQ8h and for 
many of them true masses will be measured with high pre- 
cision by ground-based and space-borne astro metric observa- 
tories , such as VLTI/PRIMA and Gaia (e.g. , lLaunhardt et al.l 
2008; ICasertano etal] 120081: ISozzettil l2009l and references 
therein). The possibility that objects with similar masses formed 
in different way and exhibit different chemical composition de- 
pending on the formation mechanism open exciting perspectives 
for such instruments. Therefore, they will represent key bench- 
marks for the calibration of the masses of substellar objects, for 
furthering our understanding of the dependence of the brown 
dwarf desert on the properties of the stellar hosts, and even for 
reaching final agreement on the actual definition of planets and 
brown dwarfs themselves based not simply on semantics but 
rather on robust physical grounds. 
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Table 3. Close companions to solar-type stars with (projected) masses between 10 My and 80 My. The first part lists the probable true 
substellar objects, the second part the objects for which only the projected masses or poor constraints on the true mass are available, 
the third part the companions with true masses above the substellar limit. In the remarks, PL refers to objects with additional 
companions of planetary mass, MULT multiple systems (additional stellar companions), HIPG stars with Hipparcos acceleration 
solutions. 
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< 94 


G0IV 


1.41 


+0.37 


Endl et al. (2004), Valenti & Fischer (2005) 


HD 156846b 


10.45 


359.51 


0.99 


0.8472 






GOV 


1.43 


+0.22 


Tamuz et al. (2008). MULT 


HD 162020b 


13.75 


8.428198 


0.072 


0.277 






K2V 


0.78 


+0.11 


Udrv et al. (2002) * 


HD 167665b 


50.3 


4385 


5.47 


0.337 






GOV 


1.11 


-0.17 


Patel et al. (2007),"valenti & Fischer (2005) 


HD 174457b 


65.8 


840.800 


1.90 


0.23 








1.07 


-0.18 


Nidever et al. (2002), Valenti & Fischer (2005) 


HD 184860b 


32.0 


693 


1.4 


0.67 






K2V 


0.77 


-0.04 


Vogt et al. (2002). Valenti & Fischer (2005). MULT 


HD191760b 


38.17 


505.65 


1.35 


0.63 






G3IV/V 


1.28 


+0.29 


Jenkins et al. (2009) 


HD202206b 


17.4 


255.87 


0.83 


0.435 




< 149 


G6V 


1.17 


+0.35 


Correia et al. (2005), Zucker & Mazeh (2001). PL 


HD283750b 


51 


1.788 




0.002 






dK5 


0.67 


+0.03 


Halbwachs et al. (2000), MULT 


TYC2534-698-lb 


39.7 


103.698 


0.44 


0.385 






G2V 


1.00 


-0.25 


Kane et al. (2009) 


HAT-P-13C 


15.2 


428.5 


1.19 


0.691 






G4V 


1.22 


+0.43 


Bakos et al. (2009). PL 


Probable stellar companions with projected mass in the substellar range 


HD18445B 


45 


554.6 




0.558 


-16.1 


187 (20) 
178 ± 20 


K2V 


0.78 


0.00 


Halbwachs et al. (2000). Valenti & Fischer (2005). MULT ' 
Zucker & Mazeh (2001) 


HD33636B 


9.3 


2128 




0.48 


4.1 ±0.1 


142+f 




1.02 


-0.13 


Bean et al. (2007) 


HD43848B 


25 


2371 


3.4 


0.69 


55 ± 12 


120«f 


GOV 


0.89 


-0.03 


this paper, HIPG 


HD110833B 


17 


271.17 




0.784 


7.76 


146(12) 
143 ± 12 


K3V 


0.72 


0.00 


Halbwachs et al. (2000). Mishenina et al. (2004) * 
Zucker & Mazeh (2001) 


HD112758B 


34 


103.258 




0.139 


10.4 


212(44) 
213 ±43 


K0V 


0.79 


-0.56 


Halbwachs et al. £2000), Mishenina et al. (2004) 
Zucker & Mazeh (2001) 


HD140913B 


43.2 


147.968 


0.55 


0.54 


16.3 


177 (73) 
181 ±74 




1.17 


+0.13 


Nidever et al. (2002), Halbwachs et al. (2000) 
Zucker & Mazeh (2001), Valenti & Fischer (2005) ' 


HD164427B 


46 


108.55 


0.46 


0.55 


8.5 


-179+106 

-"•'-85 




1.18 


+0.13 


Zucker & Mazeh (2001) 


HD169822B 


27.2 


293.1 


0.84 


0.48 


175 


< 190 
320 


G5V 


0.91 


-0.12 


Tinnev et al. (2001) 

Vogt et al. (2002), MULT^ 


HD217580B 


68 


454.66 




0.520 


25.2 


171 (14) 

172 ± 15 


K4V 


0.69 


-0.06 


Halbwachs et al. (2000), Gray et al. (2003) " 
Zucker & Mazeh (200T) 


BD-04 782B 


48 


716.68 




0.074 


12.77 


261 (15) 
257 ± 16 


K5V 


0.67 


+0.05 


Halbwachs et al. (2000). Heiter & Luck (2003) 
Zucker & Mazeh (2001) 



" Fo r spectroscopic binaries with brown dwarf candidates observed by Hipparcos and with an actual mass estimate derived by Halbwachs et al.l 
(2000), the error distribution of the secondar y masses is not Gaussian and not symmetric. We have listed here the values of M c and the accompa- 
nying one-sided errors as denned by Eq. 5 in Halbwa chs et alj d2000h . The star has both chemical abundances and kinematics of thick disk. 

b Thick disk population. 

c Transiting svstem. lTriaud et aTJ f2009) found a significant misalignment between the planetary orbital axis and the stellar rotation axis (37. 6+22 '3 
deg). 

d Transiting svstem. lwinn et all J2009^ found a significant misalignment between the planetary orbital axis and the stellar rotation axis (37.3±3.7 
deg). 

E The orbit in lPerrier et alj d2003h is preliminary, with only limited phase coverage. iPerrier et alJ d2003h also reported no detection with AO. 

' Assuming synchornous rotation and orbit aligment Fuhrmann (2008) derived = 15 and M c = 0.14*JM5 M Q for the companion. The system is 
a quadruple with another SB at a = 2.53 arcsec (21 AU) with e=0.41 from HD98230. 

g The possibility that the star is seen nearly pole-on because of its very small projected rotational velocity is widely debated in the literature. 
Halbwachs et al. (2000) derived a very un certain astrome t ric ma ss 112 Mj with an error of 103 Mj. Additional low-mass companion close to the 
substellar boundary detected at 130 AU bv lPatience et al.l d2002h . Member of the thick disk population. 

* IZucke r & Mazeh ( 2001) derived an astromeric solution that gives a companion mass of 3.0±0.94 M Q . Not compatible with the stellar properties 
of the primary. Probably spurious due to the short orbital period. 

' lHalbwachs etai] d2000h derived an astrometric mass of 182 Mj with a large error of 470 Mj. Gle bocki & Stawikowskil dl995h estimate 
i = 22 ± 10 and then M = 136 ± 39 Mj from the rotation period (synchronous with the orbit) and projected rotational velocity, assuming the 
rotation axis is perpendicular to the orbital plane. However, misaligne ments can not be ru led out, especially considering that the system has a wide 
WD companion. The system is probably associated with the Hyades ([Catalan et al. 2008). Metallicity from[Catalan et al. (20081). 

' Also resolved with direct imaging bv lBeuzit et alJ d2004). Quintuple system (see Bonavita & Desidera (2007) and references therein). 

k Purely astrometric orbital solution in Hipparcos. 



